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INTRODUCTION 

Ecological communities comprise complex networks of inter- 
acting species that influence each other via many direct and 
indirect effects. In an ecological context, direct effects are those 
that one species exerts on another species with which it directly 
interacts: e.g., plant-herbivore interactions (Agrawal et al., 2006). 
Indirect ecological effects (IEEs) are mediated via a third species 
(Wootton, 1994): e.g., host plant effects on a predator medi- 
ated via the herbivorous prey (Astles et al, 2005). Both direct 
and indirect effects are thought to play important roles in evolu- 
tionary processes via the promotion or constraint of phenotypic 
evolution (Wootton, 1994), although specific evidence for this is 
limited (Strauss and Irwin, 2004). 

Genetic diversity is the key ingredient linking ecological and 
evolutionary processes, and sufficient genetic diversity is required 
if ecological interactions (including IEEs) are to affect the evolu- 
tion of phenotypes. The discipline of community genetics seeks 
to provide a framework that allows us to investigate the influ- 
ence of intraspecific genetic variation in one, or more, species 
on the phenotype of other interacting species (Antonovics, 1992; 
Rowntree et al, 2011b). Work has shown that there is a genetic 
basis for many direct effects (e.g., Whitham et al., 2003; Johnson 
and Agrawal, 2005; Zytynska et al, 2011) and IEEs (e.g., Bailey 
et al., 2006; Johnson, 2008). We also know that variation in the 
response of a focal species to an indirect effect can have a genetic 
component (Astles et al., 2005), and that incorporation of genetic 
variation in more than one species in multi-trophic systems can 
change the direction and magnitude of an indirect effect, in terms 
of the fitness of the focal species (Tetard- Jones et al., 2007). 



When parasitic plants and aphid herbivores share a host, both direct and indirect ecological 
effects (IEEs) can influence evolutionary processes. We used a hemiparasitic plant 
(Rhinanthus minor), a grass host (Hordeum vulgare) and a cereal aphid {Sitobion avenae) 
to investigate the genetics of IEEs between the aphid and the parasitic plant, and looked 
to see how these might affect or be influenced by the genetic diversity of the host plants. 
Survival of R. minor depended on the parasite's population of origin, the genotypes of the 
aphids sharing the host and the genetic diversity in the host plant community. Hence the 
indirect effects of the aphids on the parasitic plants depended on the genetic environment 
of the system. Here, we show that genetic variation can be important in determining the 
outcome of IEEs. Therefore, IEEs have the potential to influence evolutionary processes 
and the continuity of species interactions over time. 

Keywords: indirect ecological effects, community genetics, Rhinanthus minor, Sitobion avenae, eco-evolutionary 
feedbacks 



Therefore, there is real potential for IEEs to influence evolution- 
ary processes, but, the real impact of genetic variation in multiple 
components of a trophic system, incorporating direct and indirect 
effects, remains under-explored. 

Infection by a parasitic plant can directly influence the success 
of a host plant by dramatically changing host biomass accumula- 
tion (Cameron et al., 2005). Genetically-based variation within 
host, parasitic plant or both, has the potential to change the 
outcome of the direct interaction in terms of host biomass or 
parasite size, biomass, and fecundity (Mutikainen et al., 2000; 
Ahonen et al, 2006; Rowntree et al, 201 la). Furthermore, genetic 
diversity in the host population can influence community level 
response to infection, with genetically impoverished host com- 
munities potentially being less resistant to the parasite (Spielman 
et al, 2004). Host and parasitic plants can also indirectly influ- 
ence herbivores feeding on the other plant by either increasing 
(e.g., Ewald et al., 2011) or decreasing (e.g., Marvier, 1996) fit- 
ness. Where parasitic plants and herbivores share a host, the 
indirect effects of herbivores on the parasites have not often been 
recorded, and we do not know if such indirect effects have the 
potential to shape evolutionary processes. 

We used the hemi-parasitic plant Rhinanthus minor, a grass 
host {Hordeum vulgare) and a cereal aphid (Sitobion avenae) to 
investigate the genetics of IEEs (subsequently GIEEs) from the 
aphids to the hemi-parasite. Specifically, we wanted to see if there 
was an interaction between aphid genotype and R. minor popula- 
tion that influenced the survival and reproductive success of the 
parasitic plant. Although R. minor is not a problem in agricultural 
systems due to modern management practices, it does commonly 
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infect grasses and the Rhinanthus-Hordeum experimental system 
is well-established (Seel and Jeschke, 1999; Jiang et al., 2004). We 
also looked at how levels of genetic diversity within the host com- 
munity might influence the interactions between hosts, parasites 
and herbivores. We incorporated genetically-based intra-specific 
variation in both the aphids and the parasitic plant, and used 
either genetically uniform or genetically diverse hosts. We asked: 
(1) if there are indirect effects of the aphids on the parasitic 
plants; (2) if there is a genetic basis to any indirect effects observed 
(i.e., if different genotypes of aphid differentially impact the par- 
asitic plants); (3) if there is a genetic basis for the response of an 
organism to an indirect effect (i.e., if parasite populations vary 
in response to the aphids); (4) if the presence of genetic diver- 
sity within the host plant mediates any indirect effects observed. 
Our experiment was designed to investigate whether genetic vari- 
ation in any of the component species of a tri-trophic system 
was important in defining the indirect interactions between aphid 
herbivores and parasitic plants competing indirectly for nutrients 
via a common host. It also enabled us to determine the relative 
importance of genetic variation within a species and genetic inter- 
actions among species on the survival, reproductive potential and 
biomass production of the plants. 

MATERIALS AND METHODS 

We obtained seeds of six doubled haploid barley (Hordeum vul- 
gare L.) genotypes (Morex, Steptoe, Blenheim, Kym, Oregon 
Wolfe Dominant, Oregon Wolfe Recessive) from P. Hayes 
(Oregon State University, USA), which were bulked prior to use. 
We obtained seeds from two populations of Rhinanthus minor L.: 
the first originating in Somerset from Emorsgate Seeds (Kings 
Lynn, Norfolk, UK) and the second originating in Moray near 
Inverness from Scotia seeds (Brechin, Angus, UK). We obtained 
four cereal aphid (Sitobion avenae) genotypes (DAV95, CL07, 
HF92a, HI) from Rothamstead Research, Harpenden, UK. Aphid 
genotypes were maintained and bulked on "pearl" barley to avoid 
conditioning effects. Two of the aphid genotypes were brown 
(HF92a, CL07) and two were green (HI, DAV95) and we always 
paired one brown genotype with one green genotype to facili- 
tate identification. This gave us a total of four aphid treatments 
(Hl-HF92a; H1-CL07; DAV95-HF92a; DAV95-CL07). 

We germinated barley seeds in the dark at 23° C and then trans- 
planted them into pots (15 cm dia.) filled with horticultural sand 
(Keith Singleton Horticulture, Cumbria, UK). Barley seedlings 
were planted six to a pot and positioned in a circle around the 
edge. We used two diversity treatments: a genetically uniform 
treatment, where all six barley plants in a pot were the same geno- 
type; and a genetically diverse treatment, where each barley plant 
was a different genotype. Genetically uniform pots were planted 
for each genotype and subsequent analyses combined all single 
genotype pots for the genetically uniform treatment. Treatments 
were randomly assigned to pots, as was position and genotype of 
the barley plants in the genetically diverse treatment. 

We surface-sterilized -R. minor seeds (1% v/v sodium 
hypochlorite solution, 3 min) and germinated them in the dark 
at 4°C (see Rowntree et al, 2011a for details). Rhinanthus minor 
is a generalist parasite with the propensity to attach to multiple 
hosts simultaneously (Gibson and Watkinson, 1989). Our aim 



was to start with three parasites successfully attached to one or 
more host plants. However, as not all seedlings are successful in 
attaching to a host we planted six _R. minor seedlings ( 1-2 cm radi- 
cles) per pot. Parasite seedlings were planted in a circle in the 
center and positioned approximately 2 cm away from the barley 
plants. We used two populations of R. minor and a no R. minor 
control for each barley diversity and aphid-pair treatment. Plants 
were grown in a greenhouse (15-25°C) with supplementary light- 
ing (16:8 photoperiod) and watered every day with Vi strength 
Hoagland's solution (Hoagland and Arnon, 1950). The amount 
of nutrient solution added to the pots differed over time depend- 
ing on the requirements of the plants, but was always consistent 
among pots on a single day. If the pots remained dry (i.e., if out- 
side conditions were sunny and warm), we supplemented with 
water. In total, we planted 1230 germinated R. minor seeds in 205 
pots and 1872 barley seedlings in 312 pots. 

Two weeks after planting we counted all emerged seedlings 
of R. minor. Five weeks after planting we scored R. minor plants 
for morphological characteristics indicating host attachment and 
counted the number of attached plants (see Klaren and Janssen, 
1978). At this point we reduced the number of R. minor to 
three per pot. We removed unattached plants preferentially and 
then chose which of the remaining individuals to remove ran- 
domly. Any pots without three attached _R. minor were removed 
from the study leaving a total of 277 pots, 169 with R. minor. 
This gave us 8-12 pots per Rhinanthus-aphid combination in the 
genetically diverse treatment and 8-15 pots per Rhinanthus-aphid 
combination in the genetically uniform treatments. 

Aphids were added 6 weeks after planting (±1 day). We put six 
fourth-instar or adult individuals of two aphid genotypes (one 
brown, one green) in a 6 cm petri dish (12 aphids per pot) and 
placed it open in the center of the pots. We used two aphid geno- 
types per replicate pot as one of the goals of the experiment was 
to investigate intraspecific competition among the aphids. These 
results have been previously reported (Zytynska et al, 2014) and 
are not included in our analysis here. We covered the plants and 
pots with fine mesh bags (Insectopia, Austrey, Warwickshire, UK) 
to ensure that aphids could not escape, but remained free to move 
among all plants within a treatment. Two weeks later (±1 day), 
after sufficient time to reproduce, we counted the numbers of 
each genotype of aphid on every plant in each pot. Finally, we 
removed the covers and all aphids from the plants by hand. 

Thirteen weeks after planting, we counted the number of flow- 
ering _R. minor individuals per pot to assess survival of potentially 
reproducing plants. We also counted the number of buds, flow- 
ers and seedpods (reproductive structures) per flowing plant and 
combined these as a measure of reproductive fitness. By the end 
of the experiment there were 95 pots in the R. minor treatment 
level containing a total of 206 flowering parasitic plants. Finally, 
we harvested each barley plant separately, dried them at 60° C for 
at least 2 days and measured shoot dry weight. 

DATA ANALYSIS 

We analyzed survival and fitness of individual _R. minor plants 
at a series of key life history stages: (1) between germination 
and seedling emergence (emergence data); (2) between emer- 
gence and attachment to a host (attachment data); (3) attachment 
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through to production of flowers (flowering data); and (4) the 
total number of buds, flowers and seedpods produced (Figure 1). 
The number of plants at each life history stage depended on suc- 
cessful survival of the previous stage. Therefore, the number of 
plants naturally declines for subsequent analyses. We used gener- 
alized linear mixed models with a binary distribution, a logit link 
function and pot as a random factor to analyse the survival data. 
Rhinanthus minor population (Somerset or Inverness) and barley 
diversity (uniform or diverse) were fixed factors in the analyses of 
emergence and attachment data, and interactions were included 
in the model. Fixed factors in the analysis of flowering data were 
_R. minor population, barley diversity, green aphid genotype (HI 
or DAV95) and brown aphid genotype (HF92a or CL07). The 
final total number of aphids for each pot was used as a covariate 
and all possible interactions were included in the model. We cor- 
rected for variation in the time for aphid reproduction (either 14 
or 15 days) by dividing final total number of aphids by number 
of days for each pot. The total number of buds, flowers and seed- 
pods produced by each R. minor plant that survived were analyzed 
using a single generalized linear mixed model with a negative 
binomial distribution, a log link function and pot as a random 
factor. Rhinanthus minor population, barley diversity, green aphid 
genotype and brown aphid genotype were included as fixed fac- 
tors, total number of aphids (corrected as above) was used as a 
covariate and all possible interactions were included in the model. 

Prior to analysis, individual barley shoot dry weight data 
were square-root transformed. Q-Q plots of raw data, natural 
log transformed data and square root transformed data showed 
the latter to confirm best to normality. The shoot dry weight 
from individual barley plants was analyzed with a linear mixed 
model where pot was included as a random factor. Rhinanthus 
minor presence, barley diversity, green aphid genotype and brown 
aphid genotype were included as fixed factors, total number of 
aphids (corrected as above) and the final number of R. minor 
plants per pot were included as covariates. All possible interac- 
tions were included in the model. In order to test the effect of 
barley genotype on shoot dry weight, we subdivided the data by 
diversity treatment and reanalyzed using similar statistical models 



but replacing barley diversity with barley genotype. Data were 
split and reanalyzed as the experimental design did not allow us 
to incorporate genotype into the full model. 

All analyses used the R statistical platform version 3.0.2 (R 
Core Team, 2013). Analyses of the survival and total bud, flower 
and seedpod data were performed using the glmmadmb function 
in the glmmADMB package (Fournier et al., 2012; Skaug et al., 
2013). Analyses of the barley shoot dry weight was performed 
using the lmer function in the lme4 package (Bates et al, 2014). 
Significance values were calculated from Wald chi square tests 
using the Anova function in the car package (Fox and Weisberg, 

2011) . As part of the goal of the study was to uncover any 
potential interactions among treatments, we used a full model 
approach rather than a model simplification strategy (Crawley, 

2012) . 

RESULTS 

Only 20% of germinated seeds survived to produce flowers. When 
we broke survival down by the different life stages (Figure 1), 
70% emerged above ground to produce seedlings. Emergence 
was higher in the Inverness population at 76% compared to the 
Somerset population at 64%. There was a significant effect of R. 
minor population on emergence (X 2 = 4.87, p = 0.03), but no 
effect of host diversity (X 2 = 0.008, p = 0.93) and no effect of 
the interaction between _R. minor population and host diversity 
(X 2 = 0.55, p = 0.46; see Table 1 for survival rates and Tables 
Sla and Sib for more test details). Of the plants that emerged, 
66% successfully attached to a host plant but there were no signif- 
icant effects of either _R. minor population {X 2 = 0.01, p = 0.91) 
or host diversity (X 2 = 3.23, p = 0.07) and no significant inter- 
action between the two (X 2 = 2.32, p = 0.13; see Table 2 for 
survival rates and Tables S2a and S2b for more test details). 

Aphids were added between attachment and flowering of the 
R. minor (Figure 1) and although aphids were free to move 
around all plants within a pot, we only found aphids on the barley 
and did not find any on the _R. minor plants themselves. Thus, any 
effects of the aphids on the _R. minor were occurring indirectly via 
the host plants. 

Subsequent _R. minor survival until flowering was 43% over- 
all but varied considerably across the treatment groups (Table 3). 
There was a significant four-way interaction between host diver- 
sity, -R. minor population, green aphid genotypes and brown 
aphid genotypes (X 2 = 4.20, p = 0.04; see Tables S3a and S3b for 
more test details), but no other significant effects. There were no 
significant effects of any of the factors tested on the number of 



Table 1 | Percentage of germinated R. minor seed that emerged as 
seedlings by treatment groups. 



Diversity 


R. minor population 


Number AliveiDead 


% Germination 


Uniform 


Inverness 


256:80 


76 




Somerset 


196:122 


63 


Diverse 


Inverness 


224:70 


76 




Somerset 


188:94 


67 



Values in bold are greater than the overall survival rate of 70%. 




FIGURE 1 | Key life stages of Rhinanthus minor measured: (1) between 
germination and seedling emergence; (2) between emergence and 
attachment to a host; (3) post attachment through to production of 
flowers; (4) numbers of buds, flowers, and seedpods (reproductive 
structures) produced. Conditional survival rates during each life stage are 
shown and are dependent on successful survival of the previous life stage. 
Absolute survival rates are shown in brackets. 
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Table 2 | Percentage of R. minor seedlings that successfully attached to a host plant following emergence (conditional attachment) and 
germination (absolute attachment) by treatment groups. 



Diversity R. minor population Number AliveiDead % Attachment (conditional) % Attachment (absolute) 



Uniform Inverness 


187:69 


73 




55 


Somerset 


124:72 


63 




40 


Diverse Inverness 


145:79 


65 




49 


Somerset 


122:66 


65 




44 


Values in bold are greater than the overall conditional attachment rate of 66% and the overall absolute attachment rate of 46%. 


Table 3 | Percentage of R. minor plants that survived until flowering following successful attachment (conditional survival) and germination 


1 a hcnl i i+o ciiKwii/nH H\/ troatmont nrniinc 
(dUwIULC oUI VIVGI0 U y LI 6fl LI 1 ICI 1 1 LJIUULJo. 










Diversity R. minor Green aphid 


Brown aphid 


Number Alive: 


% Survival 


/o ourvival 


population genotype 


genotype 


Dead 


(conditional) 


(absolute) 


Uniform Inverness Dav95 


CL07 


15:17 


47 


26 




HF92a 


12:24 


33 


18 


H1 


CL07 


21:24 


47 


26 




HF92a 


14:22 


39 


21 


Somerset Dav95 


CL07 


21:11 


66 


26 




HF92a 


7:16 


30 


12 


H1 


CL07 


14:10 


58 


23 




HF92a 


19:11 


63 


25 


Diverse Inverness Dav95 


CL07 


14:18 


44 


22 




HF92a 


8:21 


28 


14 


H1 


CL07 


8:25 


24 


12 




HF92a 


15:21 


42 


21 


Somerset Dav95 


CL07 


11:22 


33 


15 




HF92a 


19:11 


36 


16 


H1 


CL07 


8:14 


36 


16 




HF92a 


10:17 


37 


16 



Values in bold are greater than the overall conditional survival rate of 43% and the overall absolute survival rate of 20%. 



R. minor seedpods (see Table 4 for median values and Tables S4a 
and S4b for more test details). 

Barley shoot dry weight was significantly reduced by the pres- 
ence of _R. minor (X 2 = 17.25, p < 0.0001) and there was a 
significant negative relationship between dry weight and the 
number of R. minor plants (X 2 = 42.18, p < 0.0001). There was 
a significant positive relationship between dry weight and the 
total number of aphids in the pot (X 2 = 184.60, < 0.0001) and 
also a significant effect of brown aphid genotype where genotype 
HF92a had a greater positive impact on dry weight than genotype 
CL07 (X 2 = 4. 15, p = 0.04; see Tables S5a and S5b for more test 
details). 

Analysis of barley genotype effects was only possible by divid- 
ing the data by host plant diversity treatment. In the uniform 



diversity treatment, all of the barley plants in a single pot were 
of the same genotype and therefore in these pots it was not 
possible for other barley genotypes to influence barley shoot 
biomass. In these pots, barley genotype (X 2 = 74.07, p < 0.0001), 
brown aphid genotype (X 2 = 5.79, p = 0.02), _R. minor pres- 
ence {X 2 = 5.23, p = 0.02), the number of R. minor per pot 
(X 2 = 33.43, < 0.0001) and total aphid number (X 2 = 155.53, 
p < 0.0001) were significant factors affecting barley shoot dry 
weight and there was a significant interaction between bar- 
ley genotype, green aphid genotype and brown aphid genotype 
(X 2 = 12.49, p = 0.02). 

In the diverse barley treatment, each barley plant in a pot was 
a different genotype. Therefore, barley shoot biomass could be 
influenced by competition among barley genotypes. In these pots, 
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Table 4 | Median numbers of buds, flowers, and seedpods produced by the R. minor plants that survived until flowering by treatment groups. 



Diversity ft minor population Green aphid genotype Brown aphid genotype Median 25/75th Quartiles 



Uniform Inverness Dav95 CL07 11.0 7.0/23.0 

HF92a 8.0 4.5/14.5 



H1 CL07 13.0 7.0/20.0 

HF92a 8.0 4.25/22.5 

Somerset Dav95 CL07 15.0 12.0/19.0 

HF92a 16.0 8.0/21.5 

H1 CL07 17.5 12.25/21.75 



HF92a 8.0 6.5/14.0 



Diverse Inverness Dav95 CL07 8.5 0.5/16.25 





HF92a 


4.5 


0.75/13.0 


H1 


CL07 


14.0 


10.0/16.75 




HF92a 


12.0 


6.5/22.0 


Somerset Dav95 


CL07 


12.0 


4.0/19.0 




HF92a 


12.0 


7.5/21.0 


H1 


CL07 


10.5 


8.25/13.5 




HF92a 


3.5 


0.25/15.75 



Numbers in bold are greater than the overall median value of 12.0 (6/40). 



there was a significant effect of barley genotype (X 2 = 504.51, 
p < 0.0001), R. minor presence (X 2 = 21.54, p < 0.0001), the 
number of R. minor per pot (X 2 = 17.73, p < 0.0001) and the 
total number of aphids {X 2 = 90.65, p < 0.0001), but no effect of 
brown aphid genotype {X 2 = 0.2774, p = 0.60) on barley shoot 
dry weight. There was, however, a significant interaction between 
barley genotype and brown aphid genotype (X 2 = 22.11, p = 
0.0004; see Tables S6 and S7 for more details). 

DISCUSSION 

We collected data on the survival and reproductive fitness of 
the annual parasitic plant R. minor over four key life stages in 
the context of genetically variable multi-trophic species interac- 
tions. We found that the highest percentage of plants died post 
attachment and after addition of aphids to the host. Survival dur- 
ing this life stage was dependent on the genetic context of all 
species involved. Previously, we have demonstrated that aphid 
distribution across host plants is influenced by the presence of R. 
minor, and when present, by genetic diversity within the R. minor 
infecting the hosts (Zytynska et al, 2014). In this study we show 
that aphid genotype can have a reciprocal effect on the survival 
of the parasitic plant. Thus, our experimental system demon- 
strates that genetically-based feedbacks can occur (sensu Genung 
et al., 2011) between indirect competitors. Specifically, we show 
that the population of R. minor influences the distribution of 
aphid genotypes across the host plants, and that aphid geno- 
types, in combination with host genetic diversity, differentially 



influence the survival of individuals from two populations of 
R. minor. 

We know from previous studies that genotype by genotype 
interactions can influence the direct interactions among species 
with important consequences for the outcome of, among other 
things, plant-herbivore (Tetard-Jones et al., 2007), host-parasite 
(Lambrechts et al, 2005; Salvaudon et al, 2005; Rowntree et al., 
2011a), and plant-plant (Fridley et al., 2007) interactions. They 
also have the potential to influence the evolutionary and coevo- 
lutionary trajectories of the species involved (Thompson, 2013). 
Genetic variation within a focal species can also indirectly affect 
other species in the wider associated community (Fritz, 1995; 
Irwin, 2006; Johnson, 2008; Schadler et al., 2010) and deter- 
mine how a particular species responds to an IEE (Astles et al., 
2005). Previous studies have demonstrated indirect effects of 
herbivores on R. minor (Bass et al., 2010). However, to our 
knowledge, this is the first time that GIEEs have been shown 
to be an important factor determining the fitness of a parasitic 
plant. Moreover, in combination with our previous analysis we 
demonstrate that there are reciprocal GIEEs between the aphids 
and the parasitic plant (Zytynska et al., 2014). This means that 
GIEEs can play potentially important roles as selective agents in 
communities of interacting species (Ter Horst, 2010), and that 
there is potential for coevolution amongst indirectly interacting 
organisms. 

We found similar levels of mortality for R. minor during 
the two initial life stages we examined (germination-emergence; 
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emergence-host attachment). During these early life stages, only 
population of origin influenced R. minor survival, and then only 
during the first life stage. Population level differences do not nec- 
essarily indicate a genetic basis of an effect, because maternal and 
other non-additive effects cannot be excluded. However, when 
plants are grown in a common environment, population differ- 
ences are a strong indicator of a genetic influence on the traits 
under investigation (Krebs, 2001), in this case, early survival of _R. 
minor. 

Aphids were added to the system following attachment of the 
parasitic plant to a host and prior to flowering. Mortality was 
far greater at this stage (attachment-flowering) compared to the 
two previous life history stages, suggesting that selection on _R. 
minor was stronger post attachment to a host. This may indicate 
strong indirect competition with the aphid herbivores sharing the 
host plants, but we were unable to test for this directly, as we 
did not have treatments without aphids. Alternatively, mortality 
rates could have increased at this life stage due to the addition 
of the aphid mesh-cages to all of the pots. Previous work has 
demonstrated that R. minor is intolerant to shading and low light 
levels (Tesitel et al., 201 1), which the aphid cages will have influ- 
enced. However, although shading was a confounding factor, it 
was consistent across treatments and does not explain the differ- 
ential influence of the aphid genotypes on survival of the parasitic 
plant. Nor does it explain the decreased survival of _R. minor with 
increasing host genetic diversity. Rather, this indicates that higher 
levels of genetic diversity in the host plants likely confer host com- 
munity resistance to infection by the parasitic plant (cf. Spielman 
et al, 2004; Kaunisto and Suhonen, 2013). 

The high levels of mortality following attachment of the par- 
asite to a host plant and prior to flowering were unexpected, as 
successful host attachment should facilitate survival of the par- 
asitic plant. The most common measures of fitness of _R. minor 
in previous studies have been size, biomass or flower and seed 
production (Cameron et al, 2006; Rowntree et al., 2011a; Tesitel 
et al., 201 1), and mortality rates of individuals are rarely reported. 
We found no effect of any of the factors tested on the combined 
production of buds, flowers and seedpods and, in fact, the major- 
ity of plants died before they were able to set seed. Although seed 
set must be the ultimate measure of reproductive success for an 
individual, we suggest that as single measures, the above traits 
are rather poor indicators of fitness in _R. minor and that survival 
is also crucial to monitor. If, as is indicated here, there is dif- 
ferential survival of plant genotypes prior to flowering and seed 
set, then there is potential for GIEEs to result in diverging lev- 
els of fitness, despite the lack of difference observed in seedpod 
production. 

Host plant shoot biomass was most strongly influenced by 
the level of infection by both the parasitic plant and the aphid 
herbivores, although we also detected an effect of the brown 
aphid genotype. Previously, we have found genetically based 
variation within both _R. minor and the aphid S. avenae to influ- 
ence host plant traits (Tetard-Jones et al, 2007; Rowntree et al., 
2011a), although not in combination with each other. Much of 
our earlier work has demonstrated that the outcome of complex 
species interactions cannot be predicted by more simplistic inter- 
actions among pairs of components (Tetard-Jones et al., 2007; 



Zytynska et al, 2010) and we reach a similar conclusion here. 
We also found no evidence that the GIEEs between the aphids 
and the R. minor had any impact on the host plants themselves. 
Rather, genotype by genotype interactions between the barley 
and the aphid herbivores had a greater impact on barley shoot 
biomass. 

When biomass data from the diverse barley treatments were 
analyzed separately, we found a strong effect of barley genotype 
as well as the level of infection by the aphids and R. minor. There 
were also higher order interactions between aphid and barley 
genotypes. This suggests that in a diverse genetic population of 
hosts, interactions between individual genotypes and interacting 
species may depend on the genotypes of neighboring hosts. Mutic 
and Wolf (2007) demonstrated the importance of neighbor geno- 
type on multiple plant traits, including growth, and here we show 
that interactions among neighbor genotypes of plants have the 
potential to change the genetic interactions between a host plant 
and its aphid herbivores. 

Our experiments include multiple genotypes and popula- 
tions of multiple species, resulting in complex statistical models 
with many terms. As our intention was to test for, and identify, 
significant interactions we have not used a model simplifica- 
tion approach in our analyses (Crawley, 2012). In addition, we 
have used mixed models for hypothesis testing. This approach 
comes with the caveat that individual higher order interac- 
tions must be interpreted with caution. However, it is notable 
that our main finding that _R. minor fitness is influenced by a 
higher order interaction between R. minor population, aphids 
and host genotype mirrors the results of Zytynska et al. (2014) 
where aphids were also influenced by these same higher order 
interactions. 

In conclusion, we show that GIEEs of aphid herbivores on par- 
asitic plants have the potential to influence community structure 
by affecting plant survival rates. Our work further demonstrates, 
and supports the findings of previous studies showing the impor- 
tance of the genetic background of an ecological community 
on the performance of component species (Fridley et al, 2007; 
Whitlock et al., 2007, 2010). The value of within-species genetic 
variation on community and ecosystem processes is still largely 
unknown in systems without a single dominant long-lived plant 
(e.g., Whitham et al, 2006; Barbour et al, 2009). We show 
that variation within shorter-lived plants can influence the cur- 
rent community, and that effects on survival and reproduction 
could indeed impact the future generations and thus evolution in 
ecological communities. 
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